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tJ ■ ABSTRACT 

We study the Yaglom law, which relates the mixed third order structure 
^ , function to the average dissipation rate of turbulence, in a uniformly expanding 

*iJ i solar wind by using the two scales expansion model of magnetohydrodynamic 

lO ! (MHD) turbulence. We show that due to the expansion of the solar wind two 

. [ new terms appear in the Yaglom law. The first term is related to the decay 

Q I of the turbulent energy by nonlinear interactions, whereas the second term is 

^^ I related to the non-zero cross-correlation of the Elsasser fields. Using magnetic 

field and plasma data from WIND and Helios 2 spacecrafts, we show that at 
lower frequencies in the inertial range of MHD turbulence the new terms become 
comparable to Yaglom' s third order mixed moment, and therefore they cannot 
be neglected in the evaluation of the energy cascade rate in the solar wind. 

Subject headings: sun: solar wind - turbulence 



Introduction 



The standard theory of the solar wind (jParker Ill963l ) predicts the existence of a plasma 



flow, coming from the hot solar corona, that expands radially in the interplanetary medium 
without being further heated. In the case of a spherically symmetric adiabatic expansion. 
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the proton temperature Tp should vary as Tp r^ R ^1"^ with the distance from the sun R. 
However, observations fr om the Hehos spac ecraft from 0.3 to 1 AU showed a temperature 



radial profile Tp ^ R (iTotten et al. 1119951) . thus imp l ying; in-situ he a ting o f the solar wind 



Marsch et al.l ( 19821 ): ITu &: MarschI ( 19951) and reference s 



plasma during the expansion (see 
therein). A first possible scenario for solar wind heating was proposed by IColemanI (119681 ). 
who ascribed that phenomenon to the existence of a turbulent cascade, both in the magnetic 
and in the velocity fields, that is eventually dissipated producing internal energy. This 
has been considered an efficient mechanism at least for R < 10 AU. This idea has been 



supported by the fact that the proton temp erature is positive. 



of magnetic field and velocity fiuctuations (JBelcher fc Davis 



y cor related with the amplitude 



19711). Kqlmog orov's theory for 



stationary, isotropic, and fully developed fiuid turbulence ( JKolmogorovll 1 94 ll ) predicts for the 
one dimensional spectral energy density E{k) in the inertial interval 



E{k) = CKs'/^k-'/', 



(1) 



where Ck ~ 1-6 is the Kolmogorov constant and k is the wave number. A similar approach 
has also been extended to magnetized fiuids, as the solar w ind (jlroshnikovl Il963l : iKraichnan 



19651: 



Dobrowolnv et a 



19801 : iGoldreich fc Sridharl Il995l ) . A recent study performed by 



MacBride et al. I (120081 ) showed that the energy cascade rates derived by this method are 



not accurate enough and poorly agree with the proton heating rates. 



Assuming a model for spherical sym metric expansion of t he solar wind (jVerma et al. 



19951 ) and using proton temperature data, IVasquez et al. I (120071 ) found the following expres- 
sion for the proton heating rate per unit mass at 1 AU 



3.6 X 10~ 



%v;^ j/(Kg-s), 



(2) 



where Vsw is the solar wind velocity in km/s and the protori tempe rature Tp is measured in K. 
An alysis of the solar wind data performed by lPilipp et al. I (Il990l ) and theoretical predictions 



by iLeamon et al. I (Il999l ) suggested that in-situ electron heating rate is comparable to Ep 



Equation 
order e 



|) implies that for fast warm streams in the solar wind the heating rate is of the 
p - 10'^ J/(Kg-s), whereas for relatively slow, cool streams a typical value of the 



heating rate is Ep ~ 10^ J/ (Kg ■ s). 

An alternative method for the derivation of the energy cascade rate is to use the Yaglom 
law, which relates the third order structure function to the mean energy ca scade rate and 
represents one of the rnost fundame ntal relations in the theory of turbulence (j Yaglom 1 1 1949 



Chandrasekhar 



1967 



Frisch 



19951). Generaliza. ti on of the Yaglom law for MHD turbulence 



was de rived bvlPolitano fc Pouauet I (Il995l . Il998l ) . ISorriso-Valvo et al.l (120071 ) , ICarbone et al. 
( 20091 ) . and iMacBride et al. I ( 20081 ) used solar wind data from various satellites to derive the 



energy cascade rate by means of the Yaglom law for isotropic MHD turbulence. The latter 
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showed that the energy cascade rates obtained through th e Yaglom law are in quahtative 
agreement with the predictions of equation ([2]). Recently lOsman et al.l (120 ll[ ) used multi 
spacecraft data from the Cluster mission to study the energy cascade rate via the anisotropic 
form of the Yaglom law. It should also be noted that relatively high level of measurement 
uncertainties in the plasma data have much less influe nce on Yaglom's third ord er moment 
than on the second order moment of MHD turbulence (JGogoberidze et al. II2012I ). 



In all those previous studies homogeneity of turbulence has been assumed, while solar 
wind expansion has not been taken into account. In this Letter we study the influence of 
the solar win d expansion on the Y a glom law. Namely, using the two scale expansion model 
developed by IZhou &: Matthaeus I (Il990l ). where the 'global' variable related to the solar 
wind expansion is considered as 'slow' variable, whereas nonlinear interactions are mainly 
determined by local, 'fast', variables, we derive a modified expression for the Yaglom law. 
We show that two extra terms appear in the Yaglom law due to the solar wind expansion. 
The first extra term is related to non-WKB decay of the turbulent energy due to nonlinear 
interactions, whereas the second term is caused by the interaction between large-scale fields 
and the cross correlation of the small scale inward and outward propagating Alfven waves. 
Using data from WIND and Helios 2 spacecrafts, we show that the novel terms in the 
expansion modified Yaglom law become comparable to Yaglom's third order moment at larger 
time scales in the inertial range of solar wind turbulence and, therefore, make significant 
contribution to the estimate of the energy dissipation rate. 

The Letter is organized as follows: modified Yaglom law in the uniformly expanding 
environment is derived in Sec. 2. The solar wind data analysis is presented in Sec. 3. 
Conclusions are given in Sec. 4. 



2. Theoretical analysis 



We consider incompressible MHD turbulence in the presence of a constant magnetic field 
Bq. The Elsasser variables Z^*" = v ± h/ y/A-Kp^ the eigenfunctions of counter propagating 
Alfven wav es, are usually considered as the most fundamental variables to study MHD 
turbulence (JBiskampI |2003| ) . The dynamics of the Elsasser variables is governed by the 
incompressible MHD equations 



/^A ^ Va ■ V J Z± + (Z^ ■ V)Z± + Vp + A±V2Z+ + A^V^Z- 



0. 



(3) 



Here p is the total (hydrodynamic plus magnetic) pressure, V^ = Bo/\/47rp is the Alfven 
velocity related to the background magnetic field, p is the mass density, A^ = (z/ ± /i)/2 
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where v is the kinematic viscosity and /x is the magnetic diffusivity. Ahhough the solar wind 
plasma is weakly collisional and, consequently, turbulent fluctuations are mainly damped by 
kinetic mechanisms, in our further analysis we will include the collisional dissipation terms i n 
equation ([3]) similarly to the other studies (jPolitano &: Pouquet lll998l : ICarbone et al.ll2009l ). 
Since Yaglom's law is observed at relatively small scales in the inertial range of the solar 
wind turbulence, the precise form of the dissipation mechanism, acting at very small scales, 
seems to be unessential for the present study. 

It is well known that smooth average properties of the solar wind vary on length scales 
of the order of the heliospheric distance R. On the other hand, turbul ent fluctuations have a 
correlation length that is much smaller than R (JMatthaeus et al.ll2005l ). This scale separation 
allows for form ulation of the so-called WKB-like transpo rt equations for MHD turbulence in 
the solar wind ( iTu et al.lll984j : IZhou &: Matthaeus Ill990l ). In the framework of this approach 
the magnetic field and the velocity are considered as a sum of mean and fluctuating parts. 
The mean parts depend only on global, slowly varying coordinate (denoted by R hereafter) 
and is related to the large scale inhomogeneities in the solar wind plasma, whereas fluctuating 
(turbulent) fields depend on both the slow and a rapidly varying c oordinate, denoted here 
by r. Large and small scales can be separated, assuming ergodicity ( JMatthaeus fc Goldstein 
19821 ). by introducing an ens emble averaging operation (). Here we use the two scales 
expansion model presented in IZhou fc Matthaeus I ( Il990l ): equations for fluctuations in the 
expanding environment are 






+ (UtVa)-Vz= 



r± -^=F 



+Z^ 



vu± 



1 



1 



2 
VBo 



-V 



u 



±v. 



+ (z^ ■ V)z^ 



P 



Vp + A^V^z+ + A^V^z" 



(4) 



Here z^ denote the turbulent fluctuations of the Elsasser variables and Bq is the large scale 
mean magnetic field. This equation is derived assuming incompressibility with respect to 
the small scale variable r, i.e. Vr ■ v = Vr ■ p = 0. 

In hydrodynamics a mean or large-scale flow sweeps the small-scale eddies without 
affecting the energy transfer between length scales. In magnetohydrodynamics a mean (or 
large-scale) magnetic field Bq sweeps oppositely propagating fluctuations z~ and z+, which 
affects the energy transfer. Therefore, contrary to hydrodynamics, in case of MHD turbulence 
both linear (sweeping) and nonliear (streining) tim escales are irnportan t for th e dyna mics of 
turbulence in the inertial interval (for a review see IZhou et al. I (120041 ): IZhou I (120101 )). The 
consept of the scale separation introduced above implies that the energy cascade timescale 
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(which in general case depends on both the sweeping and straining of turbulent fluctuations) 
is much less compared to the characteristic timescale of the solar wind expansion R/U . 

For further simpliflcations of the analysis we make several additional assumptions: we 
consider a uniform radial expansion, i.e., U = ([/, 0, 0) with U = const, and a constant large 
scale mean magnetic fleld Bq; we also assume U ^ Va, which is reasonable in the case of the 
solar wind (being U ~ 400 km/s and Va ~ 40 km/s), and flnally we assume A"*" = A^ = z/. 
For a uniform expansion VU = 2U/R, so that equation (jl]) reduces to 



dz 



± 



dt 



+ U ■ Vz^ T V^ ■ Vz^ + (z^ ■ V)z^ + 



U 



P 



-Vp + z/V^(z+ + z' 



(5) 



As we see, due to the expansion, two extra terms appear in equation [5] compared to the 
standard incompressible MHD equation for the turbulent fluctuations ( JBruno fc Carbone 
20051 ) . The second term on the LHS of equation (jSI) describes convective acceleration related 
to the solar wind expansion, whereas the flfth term d escribes interaction betwee n large scale 
flelds and the cross correlation of the Elsasser fields (jZhou fc Matthaeus I Il989l ) . 



The derivation of the Yaglom law in MHD turbulence (jPolitano &: Pouauet Ill998l) ha s 
already been presented in details without the expansion effects by ICarbone et al.l ( l2009l ). 
Derivation of the Yaglom law implies the following steps: one should consider equation ([5]) 
in two different points, x' and x, thus allowing to define R = (x' + x)/2 and r = x' — x. 
Subtracting those equations in x' and x we can derive dynamical equations for the differences 
of the Elsasser fields Az='= = z^(x') — z='=(x). We then multiply the obtained equations by 
Az^, we ensemble average and take the trace. Under the assumption of local homogeneity, 
on e finally ends u p with the dynamical equation for the evolution of (|A2;^p) (equation (12) 



m 



Carbone et al.l (120091 )). It can be readily shown that the two extra terms related to the 
expansion evolve as follows: the third term on the LHS of equation ([5]) gives f/S^dA^^p), 
whereas the forth term gives U{\Az^\'^)/R — ?7(Az+ ■ Az^)/R. Therefore, the dynamical 
equation for (|Az^p) takes the form 



d_ 



\Az^\')+V-{Az^\Az^\')+^-^{R{\Az^\'))- 



U 
R 



(Az+-Az-) = 2uVmAz 



±|2 



4, 
3 



)-^V-(5±r), 

(6) 

Here e^ = i'{diZj diZ, ) (where indices indicate the summation over the vector components) 
are the pseudo energy dissipation rates of th e corresponding Elsa sser fields. Equation (El) is 
valid even for anisotropic MHD turbulence ( ICarbone et al.l 120091 ). A further simplification 
implies assumption of local isotropy. Although MHD turbulence is known to be a nisotropic 
i n the inertial range due to the influence o f the m ean magnetic fleld (see, e.g., iBiskamp 
(120031 )). study performed by IStewarz et al. 1 (120091 ) showed that the analysis of solar wind 
data based on isotropic and anisotropic versions of the Yaglom law gives very similar results 
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for the energy dissipation rate of MHD turbulence, therefore the assumption of isotropy is 
reasonable for the purposes of the presented study. In this framework the gradient and the 
Laplacian operators become 



2 d_ 
r (9r ' 



2 d_\d_ 
r dr J dr 



(7) 



Considering stationary turbulence (thus dropping the first term on the LHS of equation IQ), 
which could be a reasonable approximation in certain fast wind streams, the integration of 
equation ([6]) in the inertial range (where the influence of the dissipation term on the RHS 
can be neglected) yields 



{Az^\Az^\') + 



U 



d 



y'^{R{\Az^\'))dy 



U 



y^(Az+ ■ Az")dy 



4 ± 

--e r. 
3 



This equation represents a generalization of the Yaglom law for uniformly expanding solar 
wind. The second term on the LHS of equa tion (|8D is related to the nonlinear interactions of 
the fluctuations. Indeed, as it was shown by lMarsch fc Tul (jl990l ). in the absence of nonlinear 
interactions and cross correlations, the energy densities of counter propagating Alfven waves 
would have a WKB dependence on heliospheric distance E^[R)p[R)~^/'^{U'^ =F V^) = const 



( JMarsch &: Tul 119891 ). In the considered case of the uniform expansion {U = const, p{R) ~ 
1/i?^) and for U ^ Va, it is easy to show that for fluctuations following WKB scaling, 
namely E^{R) ~ 1/-R, the second term on the LHS in equation ([8]) vanishes. Therefore, 
this term is directly related to the decay of the turbulent energy by means of nonlinear 
interactions during radial expansion. Note that a si milar term also app ears in the modified 
Yaglom law for decaying hydrodynamic turbulence ( JDanaila et al.lll999h . which is known to 
dominate over the Yaglom's third order term for relatively large separation r. 

The second extra term related to the solar wind expansion (the third term on the 
LHS of equation ([H])) is due to the cross correlation of the Elsasser fields. The identity 
(Az^ ■ Az~) = (Af ^) — (A6^)/47rp shows that the cross correlation is proportional to the so- 
called residual energy, i.e., the difference between kinetic and magnetic energies. Taking into 
account that both solar wind data and numerical simulations of MHD turbulence show an 
excess of magnetic energy with respect to th e energy contained i n the velocity fluctuation s 
at all scales in the i nertia l range (see, e.g., ITu fc MarschI (119911 ): iBavassano et al.l ( 1l998l ): 



Gogoberidze et al. I ( 120121 )). we conclude that the third term on the LHS of equation (jS]) is 



positive. 

Finally, it is worth noting that a dimensional analysis of equation (|8]) yields for a ratio 
between the Yaglom third order term and the terms related to the expansion to be of order 
{'Vrms/L)/{U/R), where Vrms is the rms value of the velocity fluctuations. At 1 AU this 
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ratio has the order > 10 and one can conclude that expansion effects are small (or even 
negligible) in the inertial range of the solar wind turbulence. As we show below, this is not 
the case and at least for larger time scales in the inertial range these terms are of the same 
order of magnitude. The reason is that the Yaglom law contains ensemble average of third 
order quantity which does not have specific sign, and consequently, the result obtained after 
averaging is significantly reduced compared to the dimensional estimate {ts.zf\^z^\'^) ~ f^ms- 

To study the relative importance of the third order and of the solar wind expansion 
terms, we performed data analysis of the data from WIND and Helios 2 satellites. The 
decay term (second term on the LHS) contains derivative with respect to the heliospheric 
distance, so that it requires data of turbulent fluctuations at least at two different heliospheric 
distances. Due to this reason, first we use 3 sec resolution plasma and magnetic field data 
from the WIND satellite to compare the third order Yaglom and the cross correlation terms. 
Then we use 81 sec resolution data from Helios 2 collected fro r n the same corotati ng plasma 
stream at different heliospheric distances (JMarsch fc Tulll990l : ITu fc Marschlll995l ) to assess 
the possible importance of the decay term. 



3. Solar Wind Data Analysis 

In order to estimate Yaglom's third order term and the cross correlation term, we use 
magnetic field data from the Magnetic F ield Investigation (MFI) instrument on board WIND 
at 3 sec resolution (JLepping et al.lll995[ ). Density and ve locity data are provided by the three 
dimensional plasma (3DP) instrument (JLin et al.lll995l ). We use observations made during 
a quiet fast stream. The start time of the interval is 10:00 of 2008 February 02 and the end 
time is 00:00 of 2008 February 04. During this interval, the solar wind speed remained above 
550 km/s. The energy of compressive fluctuations was an order of magnitude lower than that 
of incompressible fluctuations and, consequently, magnetic and velocity fluctuations, being 
mainly Alfvenic, were dominated by the cornpone nts perpendicular to the local mean field. 
We made use of Taylor's hypothesis (JTaylor Ill938l ) r = ±t/r to relate temporal changes r in 
the observational data to spatial variations r of the turbulent fields, by using the solar wind 
speed f/ as a transformation parameter. In the last expression the plus sign corresponds 
to the case when the radial axis is directed toward the sun, for instance in Geocentric 
Solar Ecliptic (GSE) coordinates (where the x axis points toward the sun and the z axis is 
perpendicular to the plane of the earth's orbit around the sun), the minus sign corresponds 
to the case when the radial axis is directed away from the sun, as in the Radial- Tangential- 
Normal (RTN) coordinate system. Using Taylor's hypothesis, equation ([8]) can be rewritten 



as 



aY^ + D^ + M = -e'- 



(9) 



where F± = 3{Az^\Az^\'^)/4Ut is the Yaglom term, D^ = 3U J^y^^{R{\Az^\'^))dy/4Rt^ 
is the decay term, M = ~3U J^ t/^(Az+ ■ Az^)dy/4Rt^ is the cross correlation term, a is +1 
if r = Ut and —1 for r = —Ut, depending on the reference frame, as discussed above. 

The energy cascade rate of MHD turbulence e is defined as £ = (£+ + £^)/2, so that. 



aY + D + M = -e. 



(10) 



Here Y = (Y^ + Y^)/2 is the mean Yaglom term and D = {D^ + D^)/2 is the mean decay 
term. The data from the WIND satellite are provided in GSE, so that a = 1. The absolute 
value of the mean Yaglom term (solid line) and the cross-correlation term (dashed line) for the 
solar wind interval studied here are presented in Figure [H The mean energy dissipation rate 
derived from the Yagl om relation e = —Y is e fti 5000 J Kg~^ s~^ in qualitative agreemen t 
with previous studies (ISorriso-Valvo et al.ll2007l : iMacBride et al. Il2008t lOsman et al.ll201ll ). 
The mean proton heating rate for the studied interval, derived by means of equation ([2]), 
Ep = 6050 J Kg~^ s~^ is indicated by the red horizontal line in Figure 1. It has to be noted 
that the Yaglom's term is derived by averaging the third order mixed term that does not have 
a fixed sign, so that it requires much more data points for stable convergence with respect 
to the second order moment with fixed sig n. Indeed, it has bee i i showed that the st able 
convergence requires up to 10^ data points (jPodesta et al. II2009I : IStewarz et al. II2009I ). If 
this con dition is not fulfilled, the n the Yaglom law is observed only in certain solar wind 
streams (ISorriso-Valvo et al.l 120071 ). Our interval contains about 5 x 10"^ points, thus giving 
quite stable results; they are found to be in agreement with other studies. 

As it is well known the residual energy in the solar wind u sually follows power l aw sc aling 
in the inertial interval (Az"'"-Az~) ~ t"' with 7 ^ 0.7 (see, e.g.. lGogoberidze et al. 1 ( 120121 ) and 
references therein). Thus, it is expected that in the inertial range the cross- correlation term 
in equations ([8])- ([9]) should behave as M ~ t"^. On the other hand, at large scales correlations 
between fluctuations at different points are weakened, i.e., (Az^-Az") ^ 2(z"'"-z~), therefore, 
for very large time separations M — )■ const. This behaviour is clearly seen in Figure [TJ 

It can be noted in Figure [T]that at tens of minutes scales in the inertial range the Yaglom 
and the cross-correlation terms are of the same order of magnitude. The ratio of these two 
terms at the scale of 30 minutes |y|/M ?a 3.5. We studied several tens of other intervals of 
fast solar wind streams and the mean value of this ratio at 30 min scales was {\Y\/M) ^ 2.5. 
Usually, the mean energy cascade rate in the solar wind is determined as the average value of 
the dis sipation rate derived v ia the relation e = —Y for time separations from 1 minute to 2 
hours ( IMacBride et al. II2008I ). Since the cross-correlation term is always positive, neglecting 
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M can lead to an overestimation of the real energy dissipation rate by tens of percent. We 
also analysed quasi stationary intervals of slow solar wind streams and the analysis showed 
that the cross-correlation term in equation ( ITUj) is much less important for slow streams. The 
typical value of the ratio is |y|/M k. 10 — 20. These findings can explain some significant 
differences found between the energy dissipation rates obtained usin g equation (l2 l) and the 



estimation of the Yaglom law in fast solar wind streams observed by lStewarz et al. I (120091 ). 
These authors found very good agreement between the proton heating rate and energy 
cascade rate for relatively cold, slow streams of the solar wind, whereas for hot, fast streams 
the Yaglom law provided significant overestimate of the cascade rate, compared to the value 
predicted by equation ([2]). 



We also studied data from Helios 2 with 81 sec cadence coL 



ected from the same coro- 



tating stream at different heliospheric di stances (IBavassano et al.lll982l : iMarsch &: Tulll990l ). 
We used the data studied previously by IZhou fc Matthaeus I (Il990l ). The first interval is on 



1976 day 76 when the distance from the sun was 0.65 AU. The start time of the second 
interval is 00:00 on 1976 day 50 and the end time is 22:04 on 1976 day 51, when the dis- 
tance from the sun was 0.87 AU (IMarsch &: Tulll990l ). The data are provided in the RTN 
reference frame so that a = — 1 in equation ( llOp . The mean Yaglom (solid line) and the 
cross-correlation (dashed line) terms at 0.87 AU are presented in Figure |2j Their behaviour 
is similar to the one found at 1 AU in the WIND data set, although the Yaglom term in 
Figure [2] is less stable. This is due to the fact that the data set contains only about 3000 
points. The absolute value of the mean decay term D (which is actually negative due to 
decay of the turbulent energy during evolution) is given by the dash-dotted line in Figure 
|2j We used a linear approximation for the estimation of the derivative with respect to the 
heliospheric distance in the D term of equation (|9]). In particular, for any variable Q we 
assumed dR[Q{R)] ^ [Q{0.87AU) - Q {0. 65 AU)] /[0. 87 AU - 0.65AU]. As can be seen, the 
Yaglom and the decay terms are pretty much comparable within all the range of time scales 
considered. In this contest one could conclude that the use of the Yaglom law without the 
decay term cannot give a reliable estimation of the cascade rate. However, it is worth stress- 
ing that using two intervals relevant to the same corotating stream, gives, in fact, an upper 
estimate of the decay term. Indeed, although the data belong to the same corotating stream, 
the time gap between the data sets is of several weeks, so that the assumption of stationary 
turbulence is questionable (indeed if the total turbulence energy changes not only because 
of nonlinear decay but also because of non-stationarity, both of those effects would influence 
the decay term); further, the rough two points approximation of the derivative with respect 
to the heliospheric distance R can also contribute to the overestimation of the real value of 
D. 
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4. Conclusions 

In this letter we studied the Yaglom law for MHD turbulence in th e expanding solar 



wind using the two scales expansion model by IZhou &: Matthaeus I (Il990l ). We derived the 
Yaglom law modified by two novel terms, which take into account the expansion effect. 
One of them is related to the energy decay by nonlinear interactions in MHD turbulence, 
while the second one is related to the non-zero cross-correlation of the Elsasser fields. Using 
magnetic field and plasma data from WIND and Helios 2 spacecrafts, we show that for 
fast solar wind streams, at large time scales in the inertial range of solar wind turbulence, 
both the decay and the cross-correlation terms are comparable to the Yaglom's third order 
mixed moment, and, therefore, they can give a significant contribution in the assessment of 
the energy cascade rate. Thus, the disagreement between the proton heating rate estimate, 
obtained by using equation (|2]), and the Yaglom law observed for fast streams in the solar 
wind can be ascribed to the fact that the extra terms in equation (|8]) related to the expansion 
are usually neglected. 

S. P.'s research has been supported by "Borsa Post-doc FOR Calabria FSE 2007/2013 
Asse IV Capitale UmanoObiettivo Operativo M.2". S. P. and V. C acknowledge the Marie 
Curie Project FP7 PIRSES-2010-269297"Turboplasmas". 
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Fig. 1. — The absolute value of the mean Yaglom term \Y\ (solid line), the cross- correlation 
term M (dashed line) and the proton heating rate (red horizontal line) derived by means of 
equation (|2]) for the fast stream detected by the WIND spacecraft. 
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Fig. 2. — The mean Yaglom term Y (solid line), the cross-correlation M term (dashed line) 
and the proton heating rate (red horizontal line) derived by means of equation ([2]) from 
Helios 2 data at 0.87 AU. The absolute value of the decay term D at 0.65 and 0.87 AU is 
given by the dash-dotted line. See text for details. 
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